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Abstract 



We present a next-to-leading order QCD calculation of the cross section for isolated 
large-pr prompt photon production in collisions of transversely polarized protons. We 
devise a simple method of dealing with the phase space integrals in dimensional reg- 
ularization in the presence of the cos(2<I>) azimuthal-angular dependence occurring for 
transverse polarization. Our results allow to calculate the double-spin asymmetry 
for this process at next-to-leading order accuracy, which may be used at BNL-RHIC to 
measure the transversity parton distributions of the proton. 



1 Introduction 



The partonic structure of spin-1/2 targets at the leading- twist level is characterized entirely by 
the unpolarized, longitudinally polarized, and transversely polarized distribution functions /, 
Af, and Sf, respectively pQ. By virtue of the factorization theorem j2], these non-perturbative 
parton densities can be probed universally in a multitude of inelastic scattering processes, 
for which it is possible to separate ( "factorize" ) the long-distance physics relating to nucleon 
structure from a partonic short-distance scattering that is amenable to QCD perturbation 
theory. Combined experimental and theoretical efforts have led to an improved understanding 
of the spin structure of longitudinally polarized nucleons, Af, in the past years. In contrast, 
the "transversity" distributions 5f, first introduced in [Hj, remain the quantities about which 
we have the least knowledge. 

Current and future experiments are designed to further unravel the spin structure of both 
longitudinally and transversely polarized nucleons. Information will soon be gathered for the 
first time from polarized proton-proton collisions at the BNL Relativistic Heavy Ion Collider 
(RHIC) |4J. Collisions of transversely polarized protons will be studied, and the potential of 
RHIC in accessing transversity 5f in transverse double-spin asymmetries Att was recently 
examined in [3] for high transverse momentum prompt photon and jet production. Several 
other studies of A TT for these reactions have been presented in the past [13 III IE] 5 as well as 
for the Drell-Yan process (31 El EEl IH] • With the exception of the latter reaction [TQl [TT] , all 
of these calculations were performed at the lowest order (LO) approximation only. As is well 
known, next-to-leading order (NLO) QCD corrections are generally indispensable in order to 
arrive at a firmer theoretical prediction for hadronic cross sections and spin asymmetries. Only 
with their knowledge can one reliably confront theory with experimental data and achieve the 
goal of extracting information on the partonic spin structure of nucleons. 

In this paper we extend the results of 5J for isolated high-p^ prompt photon production, 
PP ~~ * lX, to the NLO of QCD. Apart from the motivation given above, also interesting 
new technical questions arise beyond the NLO in case of transverse polarization. Unlike for 
longitudinally polarized cross sections where the spin vectors are aligned with momentum, 
transverse spin vectors specify extra spacial directions, giving rise to non-trivial dependence 
of the cross section on the azimuthal angle of the observed photon. As is well-known |3j, for 
Att this dependence is always of the form cos(2$), if the z axis is defined by the direction 
of the initial protons in their center-of-mass system (c.m.s.), and the spin vectors are taken 
to point in the ±x direction. Integration over the photon's azimuthal angle is therefore not 
appropriate. On the other hand, standard techniques developed in the literature for performing 
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NLO phase-space integrations usually rely on the choice of particular reference frames that are 
related in complicated ways to the one just specified. This makes it difficult to fix $ in the 
higher order phase space integration. The problem actually becomes more severe if dimensional 
regularization techniques are used for dealing with the collinear and infrared singularities, as is 
customary. Even for the kinematically rather simple Drell-Yan process the NLO calculation for 
the cross section with transverse polarization is quite more complicated as for the unpolarized 
or longitudinally polarized cases In this paper, we will present a new general technique 
which facilitates NLO calculations with transverse polarization by conveniently projecting on 
the azimuthal dependence of the matrix elements in a covariant way. This method then allows 
us to carry out phase space integrals with standard tools known from unpolarized calculations. 

After presenting our technique and verifying that it recovers the known result for the trans- 
versely polarized NLO Drell-Yan cross section, we apply it to high-pr prompt photon produc- 
tion. We also present some first numerical calculations of the cross sections and the transverse 
spin asymmetry for this process at NLO. Here we of course have to rely on some model for the 
transversity densities, for which we make use of the Soffer inequality [T2*] . As in experiment, 
we impose an isolation cut on the photon. We find a moderate size of the NLO corrections and 
the expected reduced scale dependence of the cross section at NLO. 

2 Calculation of the NLO Corrections 
2.1 Preliminaries 

The transversity density Sf(x,fi) is defined [U El 13 UH] as the difference of probabilities for 
finding a parton of flavor / at scale fi and light-cone momentum fraction x with its spin aligned 
(f f) or anti-aligned (If) to that of the transversely polarized nucleon: 



(an arrow always denotes transverse polarization in the following). The unpolarized densities 

are recovered by taking the sum in Eq. (0). When the transverse polarization is described as 

a superposition of helicity eigenstates, 5f reveals its helicity-flip, chirally odd, nature (TJ 0. 

As a result, there is no leading- twist transversity gluon density, since helicity changes by two 

units cannot be absorbed by a spin-1/2 target [110 El 1 - The property of helicity conservation 

1 We note that a gluon density does contribute beyond leading twist where it will lead to terms in 

Att strongly suppressed by inverse powers of the photon pr- An estimate of such effects could follow the lines 
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in H5J. 
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in QCD hard scattering processes implies that there have to be two soft hadronic pieces in 
the process that each flip chirality, in order to give sensitivity to transversity. One possibility, 
which we are going to consider in the following, is to have two transversely polarized hadrons 
in the initial-state and to measure double-spin asymmetries 

\ [do-(tt) - MU)] _ d6a 



A 



TT 



[dor(TT) + <MU)] da 



(2) 



Here d5a denotes the transversely polarized cross section. Att is expected to be rather small 
for most processes jHUHUHj, since gluonic contributions are absent in the numerator while in the 
denominator they often play a dominant role 2 . Nevertheless, the LO study [5] suggests that 
the asymmetry for prompt photon production should be measurable at RHIC, provided the 
transversity densities are not too small. 

According to the factorization theorem |2j the fully differential transversely polarized single- 
inclusive cross section A+B — > 7+X for the production of a prompt photon with with transverse 
momentum p?, azimuthal angle $ with respect to the initial spin axis, and pseudorapidity rj 
reads 
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with hadron-level variables 
T „. -U 



V = l+- 



W = 



S = (P A +P B ) 2 , T=(P A -P 7 ) 2 , f/ = (P B -P 7 ) 2 , 



S S + T 

in obvious notation of the momenta, and corresponding partonic ones 



t 

v = 1 + - 
s 
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S=(p a + Pbf , t = (p a - p 7 ) 2 , U=(p h ~ p^f 



Neglecting all masses, one has the relations 

s = x a XbS , t = x a T , u = XbU 

,(<) 
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(4) 



(5) 



(6) 



The d8a^ are the LO (i = 0) and NLO (i — 1) contributions in the partonic cross sections 
for the reactions ab — > jX. fiR and fip are the renormalization and factorization scales. 



2 The only exception is the Drcll-Yan process, which however suffers from rather low rates. 
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2.2 Projection Technique for Azimuthal Dependence 



Let us consider the scattering in the hadronic c.m.s. frame, assuming both initial spin vectors 
to be in ±x direction. Then, on general grounds, for a parity-conserving theory with vector 
couplings, the ^-dependence of the cross section is constrained to be of the form cos(2$): 

= cos ( 2<l> ) ( xr^: ) • ( 7 ) 



d$ cos(2$) — — - — . (8) 



dpTdr]d& \ dprdi] l 

We may obtain (d 2 5a/dpTdr]) by integrating the cross section over $ with a cos(2$) weight: 

-__ ) = - / d$ cos(2$) 
dprdr] / tx J dpTdr/d'P 

For the lowest order contribution to prompt-photon production in Eq. (j3J) one has only the 
channel qq — > 7(7. Polarization for, say, the initial quark may be projected out by 

U(p a , S a ) tZ(p a , S a ) = -jJ a [1 + 7 5 ^ a ] , (9) 

where p a and s a are the quark's momentum and transverse spin vector, and u(p a , s a ) its Dirac 
spinor. One readily finds for the LO process 



(0) 

i t i \ -• ( (-/iff 

e , (10) 



where Cp = 4/3, iVc = 3 and e q is the fractional quark charge. 

As discussed in the Introduction, in the NLO calculation one wants to make as much use as 
possible of calculational techniques established for the unpolarized case. For a single inclusive 
cross section such as prompt photon production, the appropriate methods were developed in 
[To] . They involve integration over azimuthal angles. We therefore would like to follow a 
projection analogous to Eq. (jHJ; however, we should formulate it in a covariant way. To this 
effect, we first note that the factor cos(2$)/7r in the cross section actually results from the 
covariant expression 



•F(p 7 , S a , S b ) 



TXtU 



2 (p 7 • s a ) (p 7 • s b ) + — (s a ■ s b ) 

s 



IF 



which reduces to cos(2$)/7r in the hadronic c.m.s. frame. We may, therefore, use -F(p 7 , s a , s b ) 
instead of the explicit cos(2$)/7r. 

Even though employing J 7 ^, s a , s b ) becomes a real advantage only at NLO, let us illustrate 
its use in case of the LO cross section for the partonic reaction qq — > ^g. We there have 



where 5|M| 2 is the squared invariant matrix element for the reaction with transverse polarization 
and reads: 



5\M{q<1^19)\ 2 = {ee q gf^^- 



tu 

2 (p 7 ■ s a ) (p 7 • s b ) + — (s a ■ s b ) 



c 

2 W~™™™„U^l„;CU,f|2 



(13) 



One recognizes the factor .F(p 7 , s a , s&) emerging in 5|M | . We now multiply <5|M| by J-(p 7 , s a , s&) 
equivalent to the multiplication by cos(2$)/7r in Eq. (|8|). The resulting expression may then 
be integrated over the full azimuthal phase space without producing a vanishing result, unlike 
the case of <5|M| 2 itself. This integration may again be performed in a covariant way by noting 
first that the dependence of -F(j> 7 , s a , s&) 5|M| 2 on the spin vectors comes as (p 7 • s a ) 2 (p 7 ■ s&) 2 , 
(p-y • s a) (P7 ' s 6) ( s a ' s b) ; an d (s a -Sf,) 2 . The first two of these terms correspond to contractions with 
the tensors P 7 P 7 P 7 P 7 and p'ip", respectively. Expanding these tensors into all possible tensors 
made up of the metric tensor and the incoming partonic momenta, one finds straightforwardly 

u 



7 8s 2 



f f t 2 u 2 f 

/ dn 7 (p 7 • s a ) 2 (p 7 • s fe ) 2 = / dQ 7 — (2(s a • s b ) 2 + s 2 s 2 ) = / dO. 

y dQ 7 (p 7 • s a )(p y ■ s b )(s a ■ s b ) = - j cf0 7 ^(s a • s b ) 2 = - J d0 7 ^ , (14) 

where J df2 7 denotes integration over the photon phase space, and where we have chosen both 
spin vectors to point in the same direction. We also recall that Si -p a = Si -p b = (i = a, b) and 
s 2 = s 2 = — 1. We emphasize that after the replacements ((Hj) the whole invariant phase space 
over p 1 remains to be integrated, including the (now trivial) azimuthal part, as indicated by 
the j <if2 7 on the right hand side. This is the virtue of our method that becomes particularly 
convenient at NLO. It is crucial here that the other observed ("fixed") quantities, transverse 
momentum p T and rapidity 77, are determined entirely by scalar products (p a -p 7 ) and (p& -p 7 ). 
This allows the above tensor decomposition with tensors only made up of p a and p b and of 
course the metric tensor. 

Inserting all results, and including the azimuthal part of the dfl-y integration, we find 

AC 

(5\M( q q-> ig )\ 2 ) = (ee q g) 2 —^, (15) 

i\c 

and hence, using Eq. (fT2|) . we recover Eq. (jTUJl . 

In the NLO calculation, one has 2^3 reactions ab — > ^cd. For an inclusive photon 
spectrum, one integrates over the full phase spaces dfl c and dQ d of particles c and d, respectively. 
The momentum of particle d may be fixed by momentum conservation, and the integration is 
trivial. One then ends up with 

J dfi 7 j dn c F(p^s a ,s b )5\M(ab -> jcd)\ 2 . (16) 



Besides scalar products of the s 4 (i = a, b) with p 7 , the integrand may contain terms oc (s a ■ 
Pc){sb ■ Pc) and oc (si ■ p c ). As before, we may expand the ensuing tensor and vector integrals 
in terms of the available tensors. As far as the integration over dQ c is concerned, such tensors 
may be made up of the metric tensor, p a , p b , and p 1 . It is also important to keep in mind that 
in the NLO calculation we will need to use dimensional regularization due to the presence of 
singularities in the phase space integrations. We find in d = 4 — 2e dimensions: 
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(19) 



After scalar products involving p c with the have been eliminated in this way, only those with 
(p 7 • remain. As in our LO example, when we apply the factor jF(p 7 , s a , s^), these terms 
enter as (p 7 • s a ) 2 (p 7 • s b ) 2 and (p 7 • s a )(p 7 • s b ). We then may use Eq. (fTl^) after appropriate 
modification to d = 4 — 2e dimensions: 

dfly (p 7 • s a ) 2 (p 1 ■ s b ) 2 = 
dn 7 (p 7 - s a )(p 1 - s b )(s a - s b ) = - I dVt 



da 



t 2 u 2 



7 A(l-e)(2-e)s 2 
tu 
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s b y 



(20) 



2(1 -e)s y 

After this step, there are no scalar products involving the S{ left in the squared matrix element 
(except the trivial s a ■ s b = —1). We may now integrate over all phase space, employing 
techniques familiar from the corresponding calculations in the unpolarized and longitudinally 
polarized cases. As a check, we have applied our method to the Drell-Yan transversity cross 
section and recovered the known NLO result ^T] in a straightforward manner. For the interested 
reader, we list some details of this calculation in the Appendix. 



2.3 Details of the NLO Calculation for Prompt Photon Production 

From here on, all steps in the calculation are fairly standard, albeit still involved and lengthy. 
Since many of them have been documented in previous papers [IBJ H3 HH1 HH1, we only give a 
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brief summary here. We emphasize that the general method we have employed is to perform 
the integrations over the phase space of the unobserved particles in the 2^3 contributions 
analytically. We have also simultaneously calculated the unpolarized cross section and found 
agreement with the expressions available in the literature |17| ITH] . 

At NLO, there are two subprocesses that contribute for transverse polarization: 

qq -> , 

qq - 7 X. (21) 

The first one of course was already present at LO, where X = g. At NLO, one has virtual 
corrections to the Born cross section (X = g), but also 2 — * 3 real emission diagrams, with 
X = gg+qq+q'q' ■ For the second subprocess, X = qq. All contributions are treated as discussed 
in the previous subsection, i.e., we project on their cos(2<3>) dependence by multiplying with the 
function .F(jo 7 , s a , s&) in Eq. (JTTJ) and integrating over the azimuthal phase space using Eqs. (fTTjl 
and (JUJ). 

Owing to the presence of ultraviolet, infrared, and collinear singularities at intermediate 
stages of the calculation, it is necessary to introduce a regularization. Our choice is dimensional 
regularization, that is, the calculation is performed in d = 4 — 2e space-time dimensions. 
Subtractions of singularities are made in the MS scheme throughout. 

Projection on a definite polarization state for the initial partons involves the Dirac matrix 75, 
as is evident from Eq. Q. It is well known that dimensional regularization becomes a somewhat 
subtle issue if 75 enters the calculation, the reason being that 75 is a genuinely four- dimensional 
object with no natural extension to d 7^ 4 dimensions. Extending the relation {75, 7 M } = 
to d dimensions leads to algebraic inconsistencies in Dirac traces with an odd number of 75 
[2*Uj . Owing to the chirally odd nature of transversity, in our calculation all Dirac traces 
contain two 75 matrices, and there should be no problem using a naive, totally anticommuting 
75 in d dimensions. Nevertheless, we also did the calculation using the widely-used "HVBM 
scheme" j^I] for 75, which is known to be fully consistent. It is mainly characterized by splitting 
the c?-dimensional metric tensor into a four- dimensional and a (d — 4)-dimensional one. In 
the four-dimensional subspace, 75 continues to anti-commute with the other Dirac matrices; 
however, it commutes with them in the (d — 4)-dimensional one. The HVBM scheme thus leads 
to a higher complexity of the algebra 3 and of phase space integrals. We found the same final 
answers for both 75 prescriptions in all our calculations. 

Ultraviolet poles in the virtual diagrams are removed by the renormalization of the strong 
3 We use the program Tracer to perform Dirac traces in d dimensions. 
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coupling constant at a scale /xr. Infrared singularities cancel in the sum between virtual and 
real-emission diagrams. After this cancellation, only collinear poles are left. These result for 
example from a parton in the initial state splitting collinearly into a pair of partons, corre- 
sponding to a long-distance contribution in the partonic cross section. From the factorization 
theorem it follows that such contributions need to be factored, at a factorization scale /ip, into 
the parton distribution functions. A similar situation occurs in the final-state. The high-pr 
photon may result from collinear radiation off a quark, which again is singular. This singular- 
ity is absorbed into a "quark-to-photon" fragmentation function [TT1 ITS] that describes photon 
production in jet fragmentation and hence by itself contains long-distance information. The 
fragmentation contribution has not been written down in Eq. It has a structure similar to 
Eq. (JHJ), but with an extra integration over the fragmentation function. Its size also depends 
on the experimental selection of prompt photon events, as we will discuss below. 

The subtraction of initial-state collinear singularities is particularly simple in case of transver- 
sity since there is no gluon transversity and only q — > qg collinear splittings can occur. Only 
the process qq — > ■ygg has such poles. Their cancellation is effected by adding a "counterterm" 
that has the structure (for radiation off the initial quark) 

f 1 dx 5Hqq (x, M dSa ^^ s ' xt > M ' £ ) S (x (s + t)+ u) , (22) 

where in the MS scheme 



5H qg (z,fi F ) =(-l+ jE - ln47rj 5P qq {z) (j^- ) , (23) 



with the LO transversity splitting function |2Sj 



SP„(z) = C, 



1 



Here the "plus" -distribution is defined in the usual way. As indicated in Eq. (J2"2"|) . the 2-^2 
cross section in the integrand needs to be evaluated in d dimensions. The result, which turns 
out to be the same in the anticommuting 75 and the HVBM schemes, is given by 

' d5a { q °^ g (s,t,u,e) \ _ 2C F aa s 2 /i 2e f Air fi 2 s Y 2{1 - e + e 2 ) f e 2 s 2 \ 
It / ~ "iVc ~^ 6q r(l-e) \ (1 -e)(2-e) \ ~ £ ~ 2tl ) ' 

(25) 

Needless to say that we have applied also here our "projector" .F(p 7 , s a , s&) of Eq. (jTTj) and 
performed the integration over the scalar products involving spin vectors according to Eq. (120(1 . 

In the final-state collinear case, an expression very similar to Eq. (|22|) is to be used, involving 
now the unpolarized quark-to-photon splitting function 

P iq {z) = 1 + (1 " Z? (26) 



and the 2 — > 2 "pure-QCD" transversity cross sections in d dimensions, given by: 
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(27) 



In these expressions, we have neglected contributions oc 0(e 2 ), which do not contribute. Then, 
the results for a fully anticommuting 75 and for the HVBM prescription are again the same. 

Before coming to our final results, we would like to make two more comments on the use of 
our "projector" on the azimuthal-angular dependence, Eq. ([lip. In an NLO calculation, carried 
out in d dimensions, we could have a projector that by itself contains terms oc e. Indeed, some 
of the Born cross sections, when evaluated in d dimensions, suggest a projector of the form 



Sb 



TttlL 



tu 

2 (p 7 • s a ) (p 7 • s b ) + (1 - ae)— (s a ■ s b ) 

s 



:28) 



with some constant a. Clearly, the final answer of the calculation must not depend on a because 
our projection is a physical operation which could be done in experiment. We have used the 
above projector with an arbitrary a and checked that indeed no answer depends on a. Also, we 
have integrated all squared matrix elements over the spin vectors without using any projector 
at all. This amounts to integrating cos(2$) over all < $ < 2tt, and, as expected, we get zero 
in the final answer. It should be stressed, however, that individual pieces in the calculation 
(the virtual, the 2 — > 3, and the factorization part) do not by themselves integrate to zero, 
but only their sum does. In this way, we have a very powerful check on the correctness of our 
calculation. 

2.4 Final Results for Inclusive and Isolated Photon Cross Sections 



For both subprocesses, the final results for the NLO corrections can be cast into the following 
form: 



d5a { ^ jX (s,v,w,n R ,n F ) \ _ aa s {fi R ) 



dvdw 
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+ E\nw + Flnv + Gln(l - v) + H ln(l - w) + Jln(l - vw) 



+Jln(l -v + vw) + K 
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ln(l — v + vw) 



1 — w 



(29) 



where all coefficients are functions of v and w, except those multiplying the distributions 5(1 — 
w), 1/(1 — w) + , [ln(l — w)/(l — w)] + which may be written as functions just of v. Terms with 
distributions are present only for the subprocess qq — > 7X. The coefficients in Eq. (129)1 are too 
lengthy to be given here but are available upon request. 

Let us now specify the function X &nal (z = 1 — v + vw). It results from the configurations 
where the photon is collinear with a final-state quark or antiquark. As we discussed earlier, these 
will lead to final-state collinear singularities that are absorbed, at the factorization scale 4 /i^, 
into photon fragmentation functions. The actual form of Z final depends on the kind of photon 
signal under consideration. Let us first consider the fully inclusive cross section. In this case, 
one just counts all photon candidates in the kinematical bin, without imposing any constraint 
on additional particles in the event. This is the simplest cross section and the one usually 
measured in fixed-target experiments. In the theoretical calculation, final-state singularities 
arise and there is a need to introduce a fragmentation contribution, as discussed earlier. 

At collider energies, the background from pions decaying into photon pairs is so severe 
that so-called isolation cuts are imposed on the photon. The basic idea is that photons that 
have little hadronic energy around them are less likely to result from ir° decay. The standard 
procedure is to define a "cone" around the photon by a/ (A77) 2 + (A0) 2 < R, where typically 
R « 0.4 . . . 0.7, and to demand that the hadronic transverse energy in the cone be smaller than 
rpr, where r is a parameter of order 0.1. For the theoretical calculation, isolation implies a 
strong reduction of the size of the fragmentation contribution because photons produced by 
fragmentation are always accompanied by a certain amount of hadronic energy. A slightly 
refined type of isolation has been proposed in j21]. Again a cone is defined, centered on the 
photon, within which the hadronic transverse energy must not exceed the limit r px- However, 
one chooses a larger r ~ 1 and then further restricts the hadronic energy by demanding that for 
any r < R the hadronic energy inside a cone of opening r be smaller than roughly r(r / R) 2 Pt- 
In other words, the closer hadronic energy is deposited to the photon, the smaller it has to be 
in order for the event to pass the isolation cut. This isolation method has not yet been used 
in any experiment, but it is possible that it will become the choice for the Phenix experiment 
at RHIC [23]. On the theoretical side, it has the advantage that it "eliminates" any kind of 
4 We could also choose a final-state factorization scale \j! f ^ hf here. 
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fragmentation contribution j21] because fragmentation is assumed to be a (mainly) collinear 
process, and no hadronic activity is allowed exactly parallel to the photon. 



We recall from the previous section that we have performed an analytical integration over 
the full phase space of the unobserved particles in the final-state. This seems at first sight 
to preclude the implementation of an isolation cut "afterwards". However, as was shown in 
[211 1211 5 ^ is possible to impose the isolation cut in an approximate, but accurate, analytical way 
by introducing certain "subtraction cross sections". The approximation is based on assuming 
the isolation cone to be rather narrow. In this case, dependence on the cone opening can be 
shown to be of the form a ln(R) +b + 0(R 2 ). a and b are straightforwardly determined and yield 
a very accurate description of isolation even at R = 0.7. Analytical calculations [213 12Z| are 
therefore as capable to describe the isolated prompt-photon cross section as NLO computations 
in which phase space integrals are performed numerically employing Monte-Carlo techniques 

EH El EH. 



For the cases of the fully-inclusive ("incl.") cross section, the standard isolation ("std."), 
and for the isolation proposed in [21] ("smooth") the function Z final (z = 1 — v + vw) takes the 
following forms: 



' P^(z) ]n(&) incl. 



P iq {z) In (^J + 0(1 - z[l + t]) [P jq (z) In [ {1 ~ z) *f T * ) + z\ std. (30) 
P iq {z) In ( (1 ~f T f R2 ) smooth . 



One can see the presence of the quark-to-photon splitting function P 7g of Eq. (|26|). as is expected 
for contributions resulting from near-collinear photon emission in the final-state. It also becomes 
clear that for the standard isolation the dependence on the final-state factorization scale is 
reduced and disappears altogether for the isolation of |23]. This is in line with our remarks 
above about the size of the fragmentation contribution in these cases. 



3 Numerical Results 



In this Section, we present a first numerical application of our analytical results. We focus on the 
main features of the NLO corrections and describe their impact on the cross section ddcr/dp? 
and the spin asymmetry A^ T . Our predictions will apply for prompt photon measurements 
with the Phenix detector at RHIC. This implies that the pseudorapidity region \rj\ < 0.35 is 
covered, and only half of the photon's azimuthal angle. Using Eq. (jZJ) we restore the cos(2$) 
dependence of the cross section. We take the two quadrants in $ covered by the Phenix 
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detector to be —it /A < $ < 7r/4 and 3tt/A < $ < 5tt/4 and integrate over these. This gives 
( + Is^/f) cos(2$)<i$ = 2. We consider photons isolated according to the isolation of [21] 
discussed above, using R = 0.4 and r = 1. 

Before we can perform numerical studies of A^ T we have to model the 5f we will use. 
Nothing is known experimentally about transversity so far. The only guidance is provided by 
the Soffer inequality [TJ] 

2\Sq(x)\<q(x)+Aq(x) (31) 

which gives an upper bound for each Sf. As in jH] we utilize this inequality by saturating the 
bound at some low input scale /i — 0.6 GeV using the NLO (LO) GRV [21] and GRSV ("stan- 
dard scenario") |3U] densities q(x,fj,o) and Aq(x,fj,o), respectively. For \i > fi the transversity 
densities 5f(x,fj,) are then obtained by solving the evolution equations with the LO 0123] or 
NLO [HI EI] kernels. Obviously, the sign to be used when saturating the inequality is at our 
disposal; we choose all signs to be positive. We refer the reader to [5J for more details on our 
model distributions. We note that we will always perform the NLO (LO) calculations using 
NLO (LO) parton distribution functions and the two-loop (one-loop) expression for a s . 

Figure ^ shows our results for the transversely polarized prompt photon production cross 
sections at NLO and LO for two different c.m.s. energies. The lower part of the figure displays 
the so called "if-factor" 

K = lt5a^- (32) 



One can see that NLO corrections are somewhat smaller for y S — 500 GeV and increase with 
Pt- As we have mentioned in the Introduction, one reason why it is generally important to know 
NLO corrections is that they should considerably reduce the dependence of the cross sections 
on the unphysical factorization and renormalization scales. In this sense, the i^-factor has 
actually limited significance since it is likely to be rather scale dependent through the presence 
of the LO cross section in its denominator. The improvement in scale dependence when going 
from LO to NLO is, therefore, a better measure of the impact of the NLO corrections. The 
shaded bands in the upper panel of Fig. ^ indicate the uncertainties from varying the scales 
in the range pr/2 < = — ^Pt- The solid and dashed lines are always for the choice 
where all scales are set to pt, and so is the K factor underneath. One can see that the scale 
dependence indeed becomes much weaker at NLO. 

Figure [2] shows the spin asymmetry which is perhaps the main quantity of interest 

here, calculated at LO 5 (dashed lines) and NLO (solid lines). We have again chosen all scales 

5 We note that our LO asymmetries are larger than those reported in [5]. This is due to an error in the 
numerical computation in [5]. Our LO curves in Fig. 2 correct this mistake. 
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Q ^ L_l l I l l l l I l l l l I l l 

10 20 30 

P T [GeV] 

Figure 1: Predictions for the transversely polarized prompt photon production cross sections at 
LO and NLO, for y/S = 200 and 500 GeV. The LO results have been scaled by a factor of 
0.01. The shaded bands represent the theoretical uncertainty if /xp (= (xr) is varied in the range 
Pt/2 < ix f < 2p T . The lower panel shows the ratios of the NLO and LO results for both c.m.s. 
energies. 



to be pt- Due to a larger K factor for the unpolarized cross section, the asymmetry is smaller 
at NLO than at LO. We also display in Fig. |2]the statistical errors expected in experiment. 
They may be estimated by the formula 0] 

5ArL T ~ - 1 , (33) 

where P is the transverse polarization of each beam, £, the integrated luminosity of the collisions, 
and o"bin the unpolarized cross section integrated over the pr-bin for which the error is to be 
determined. We have used P = 0.7 and C = 320 (800) /pb for VS = 200(500) GeV. 
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Figure 2: Predictions for the transverse spin asymmetry A^ T for isolated prompt photon production 
in LO and NLO for y/S = 200 and 500 GeV. The "error bars" indicate the expected statistical 
accuracy for bins in p T (see text). 



4 Conclusions 



We have presented in this paper the complete NLO QCD corrections for the partonic hard- 
scattering cross sections relevant for the spin asymmetry A^ T for high-p^ prompt photon pro- 
duction in transversely polarized proton-proton collisions. This asymmetry could be a tool to 
determine the transversity content of the nucleon at RHIC. 

Our calculation is based on a largely analytical evaluation of the NLO partonic cross sections. 
We have presented a simple technique for treating, in an NLO calculation, the azimuthal-angle 
dependence introduced by the transverse spin vectors. We will apply this technique to other 
Att in the future, such as for inclusive pion and jet production 

We found that at RHIC energies the NLO corrections to the polarized cross section are 
somewhat smaller than those in the unpolarized case. The transversely polarized cross section 
shows a significant reduction of scale dependence when going from LO to NLO. 
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Appendix: NLO Transversity Drell-Yan Cross Section 
with Projection Technique 



In this Appendix we briefly report the results we find for the NLO corrections to the Drell-Yan 
"coefficient function" SC DY when using our projection method of Sec. 2.2. For details on the 
kinematics for the process, see [TU1 ITT] . We use a fully anticommuting 75 and choose the scales 
A*f — f^R — Q everywhere, with Q the dilepton mass. The LO cross section and the virtual 
corrections at NLO rely on the underlying 2^2 reaction qq —>■ l + l~ . The real-emission NLO 
2 — > 3 process is qq — > l + l~g. We apply our projector, Eq. (JTTJ), to the squared matrix elements 
for each of these processes and integrate over the appropriate phase spaces. For the 2^3 
process this gives: 



^3 



,2e 



2nT(l - 2e) 
'ln(l - 



+8z 



13 _ 

3e 3 
In ,2 

Az- 



1 - z 



29 
18 

6z- 



, 4 26 
5(l-z)+[ — -- 



(1 



\n 2 z 



where z = Q 2 / s. For the virtual contributions we get 



6C: 



DY 



a s Cf(4:7t 



\2e 



2nT(l -2e) 

and for the MS collinear-factorization term 

a s C F (47r) 2£ 



22 



3s 



2ttT(1 - 2e) LV^ 2 



2 + Hu(i_ z ) 



+ 4(1 



116 
4 

- + 

e 



8(1 



26 

y 



l-z) 



(34) 



(35) 



(36) 



Adding all terms, the poles cancel, and one obtains the NLO MS coefficient function: 



SC BY (z) 



2tt 



-71 

3 



8 5(1 - z) 



ln(l - z) 
1 - z 



a ln z 
— Az bz 

l-z 



\n 2 z 
l-z 



+ 4(1 



(37) 



15 



in agreement with |TT] . 



References 

[1] R.L. Jaffe and X. Ji, Phys. Rev. Lett. 67, 552 (1991); Nucl. Phys. B375, 527 (1992). 

[2] S.B. Libby and G. Sterman, Phys. Rev. D18, 3252 (1978); 

R.K. Ellis, H. Georgi, M. Machacek, H.D. Politzer, and G.G. Ross, Phys. Lett. 78B, 281 
(1978); Nucl. Phys. B152, 285 (1979); 

D. Amati, R. Petronzio, and G. Veneziano, Nucl. Phys. B140, 54 (1980); Nucl. Phys. 
B146, 29 (1978); 

G. Curci, W. Furmanski, and R. Petronzio, Nucl. Phys. B175, 27 (1980); 

J.C. Collins, D.E. Soper, and G. Sterman, Phys. Lett. B134, 263 (1984); Nucl. Phys. 

B261, 104 (1985); 

J.C. Collins, Nucl. Phys. B394, 169 (1993). 
[3] J. P. Ralston and D.E. Soper, Nucl. Phys. B152, 109 (1979). 

[4] See, for example: G. Bunce, N. Saito, J. Soffer, and W. Vogelsang, Annu. Rev. Nucl. Part. 
Sci. 50, 525 (2000). 

[5] J. Soffer, M. Stratmann, and W. Vogelsang, Phys. Rev. D65, 114024 (2002). 

[6] K. Hidaka, E. Monsay, and D. Sivers, Phys. Rev. D19, 1503 (1979); 
X. Ji, Phys. Lett. B284, 137 (1992). 

[7] X. Artru and M. Mekhfi, Z. Phys. C45, 669 (1990). 

[8] R.L. Jaffe and N. Saito, Phys. Lett. B382, 165 (1996). 

[9] J.L. Cortes, B. Pire, and J.P Ralston, Z. Phys. C55, 409 (1992); 
C. Bourrely and J. Soffer, Nucl. Phys. B445, 341 (1995); 
V. Barone, T. Calarco, and A. Drago, Phys. Rev. D56, 527 (1997); 
O. Martin and A. Schafer, Z. Phys. A358, 429 (1997); 
M. Miyama, Nucl. Phys. (Proc. Suppl.) 79, 620 (1999). 

[10] W. Vogelsang and A. Weber, Phys. Rev. D48, 2073 (1993); 

A. P. Contogouris, B. Kamal, and Z. Merebashvili, Phys. Lett. B337, 169 (1994); 

B. Kamal, Phys. Rev. D53, 1142 (1996); |hep 1 ph/9807217l (unpublished); 



16 



O. Martin, A. Schafer, M. Stratmann, and W. Vogelsang, Phys. Rev. D57, 3090 (1998); 
D60, 117502 (1999). 

[11] W. Vogelsang, Phys. Rev. D57, 1886 (1998). 

[12] J. Soffer, Phys. Rev. Lett. 74, 1292 (1995); 
D. Sivers, Phys. Rev. D51, 4880 (1995). 

[13] A comprehensive review on transversity can be found in: V. Barone, A. Drago, and P.G. 
Ratcliffe, Phys. Rept. 359, 1 (2002). 

[14] X. Ji, Phys. Lett. B289, 137 (1992). 

[15] J. Soffer and O.V. Teryaev, Phys. Rev. D56, 1353 (1997). 

[16] R.K. Ellis, M.A. Furman, H.E. Haber, and I. Hinchliffe, Nucl. Phys. B173, 397 (1980); 
D.W. Duke and J.F. Owens, Phys. Rev. D26, 1600 (1982); E: D28, 1227 (1983). 

[17] P. Aurenche, A. Douiri, R. Baier, M. Fontannaz, and D. Schiff, Phys. Lett. B140, 87 
(1984); 

P. Aurenche, R. Baier, M. Fontannaz, and D. Schiff, Nucl. Phys. B297, 661 (1988). 
[18] L. E. Gordon and W. Vogelsang, Phys. Rev. D48, 3136 (1993). 

[19] B. Jager, A. Schafer, M. Stratmann, and W. Vogelsang, Phys. Rev. D67, 054005 (2003). 

[20] M. S. Chanowitz, M. Furman, and I. Hinchliffe, Nucl. Phys. B159, 225 (1979). 

[21] G. 't Hooft and M. Veltman, Nucl. Phys. B44, 189 (1972); 

P. Breitenlohner and D. Maison, Commun. Math. Phys. 52, 11 (1977). 

[22] M. Jamin and M.E. Lautenbacher, Comput. Phys. Commun. 74, 265 (1993). 

[23] F. Baldracchini et al, Fortschr. Phys. 30, 505 (1981); 
J. Blumlein, Eur. Phys. J. C20, 683 (2001); 

A. Mukherjee and D. Chakrabarti, Phys. Lett. B506, 283 (2001). 
[24] S. Frixione, Phys. Lett. B429, 369 (1998). 

[25] A. Bazilevsky, Phenix Collab., Proceedings of the "Circum- Pan-Pacific RIKEN Sympo- 
sium High Energy Spin Physics", Nov. 3-6 1999, Wako, Japan, RIKEN Review 28, 15 
(2000). 



17 



[26] P. Aurenche, R. Baier, and M. Fontannaz, Phys. Rev. D42, 1440 (1990); 

E. L. Berger and J.-W. Qiu, Phys. Lett. B248, 371 (1990); Phys. Rev. D44, 2002 (1991); 
L. E. Gordon and W. Vogelsang, Phys. Rev. D50, 1901 (1994). 

[27] S. Catani, M. Fontannaz, J.-P. Guillet, and E. Pilon, JHEP 0205, 028 (2002). 

[28] H. Baer, J. Ohnemus, and J. F. Owens, Phys. Rev. D42, 61 (1990). 

[29] M. Gluck, E. Reya, and A. Vogt, Eur. Phys. J. C5, 461 (1998). 

[30] M. Gluck, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D63, 094005 (2001). 

[31] S. Kumano and M. Miyama, Phys. Rev. D56, 2504 (1997); 

A. Hayashigaki, Y. Kanazawa, and Y. Koike, ibid. D56, 7350 (1997). 

[32] A. Mukherjee, M. Stratmann, and W. Vogelsang, work in progress. 



18 



